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Abstract: We identify and investigate novel multi-lepton signatures of extended Higgs sec- 
tors at the LHC in the guise of CP- and flavor-conserving two-Higgs-doublet models (2HDMs). 
Rather than designing individual searches tailored to specific 2HDM signals, we employ the 
combination of many exclusive multi-lepton search channels to probe the collective signal 
from the totality of production and decay processes. Multi-lepton signals of 2HDMs can arise 
from a variety of sources, including Standard Model- like production of the CP-even scalars, h 
and H, through gluon- fusion with h,H — >■ ZZ^*\ or associated production with vector bosons 
or top quarks, with h,H -> WW ( *\ZZ ( *\tt. Additional sources include gluon-fusion pro- 
duction of the heavy CP-even scalar with cascade decays through the light CP-even scalar, 
the CP-odd scalar, A, or the charged scalar, H^, such as H — > hh, H — > AA, H — > H + H~, 
H ZA, with A Zh, tt, H± Wh, and h ->• WW*,ZZ*,tt. Altogether, the combined 
multi-lepton signal may greatly exceed that of the Standard Model Higgs boson and provides 
a sensitive probe of extended Higgs sectors over a wide range of parameters. As a proof of 
principle, we use a factorized mapping procedure between model parameters and signatures 
to determine multi-lepton sensitivities in four different flavor conserving 2HDM parameter 
spaces by simulating the acceptance times efficiency in 20 exclusive multi-lepton channels 
for 222 independent production and decay topologies that arise for four benchmark 2HDM 
spectra within each parameter space. A comparison of these sensitivities with the results of a 
multi-lepton search conducted by the CMS collaboration using 5 fb^ 1 of data collected from 
7 TeV pp collisions yields new limits in some regions of 2HDM parameter space that have not 
previously been covered by other types of direct experimental searches. 
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1 Introduction 

Probing the mechanism of electroweak symmetry breaking (EWSB) is one of the primary 
objectives of the Large Hadron Collider (LHC). Fulfilling this goal includes characterization 
of the Standard Model-like Higgs boson corresponding to excitation of the scalar condensate 
responsible for EWSB [1, 2]. Yet it also extends much more broadly to include the search for 
additional Higgs states that could be a window into the underlying physics of EWSB. 

Two Higgs doublet models (2HDMs) offer a canonical framework for extended electroweak 
symmetry breaking. Indeed, in many extensions of the minimal Standard Model (SM), super- 
symmetric or otherwise, the Higgs sector is extended to two scalar doublets [3]. It is therefore 
worthwhile to study the generic features of the 2HDM scenario independent of the specific 
underlying model, purely as an effective theory for extended EWSB. The phenomenology of 
2HDMs is rich, as five physical Higgs sector particles remain after EWSB: two neutral CP- 
even scalars, h, H; one neutral CP-odd pseudoscalar, A; and two charged scalars, H + and 
H~ . All of these states could have masses at or below the TeV scale, in a regime accessible to 
the LHC. The parameter space of the 2HDM scenario is large enough to accommodate a wide 
diversity of modifications to the production and decay modes of the lightest Higgs boson, 
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as well as to provide non-negligible production mechanisms for the heavier Higgs states that 
may decay directly to SM final states, or through cascades that yield multiple Higgs states. 

Much of the study of 2HDM phenomenology to date has been devoted to the specific 
setup that arises in minimal supersymmetric models [4], which occupies a restricted subset 
of possible 2HDM signals. Even more general 2HDM studies [5-7] have largely focused on 
the direct production and decays of scalars in SM-like channels, or on specific cascade decays 
between scalars. In this work, we wish to pursue a more inclusive objective: the sensitivity 
of the LHC to the sum total of production and decay modes available in a given 2HDM, 
including both direct decays of scalars and all kinematically available scalar cascades. Such 
an approach exploits the large multiplicity of signals arising from production and decay of 
the various states in an extended EWSB sector. 

Searches for final states involving three or more leptons are well matched to this objec- 
tive, since both direct scalar decays and scalar cascades populate multi-lepton final states 
with low Standard Model backgrounds. The CMS multi-lepton search strategy [8, 9] is par- 
ticularly well-suited in this respect, since its power lies in the combination of numerous ex- 
clusive channels. While the sensitivity to new physics in any individual channel alone is not 
necessarily significant, the exclusive combination across multiple channels can provide con- 
siderable sensitivity. This is particularly effective in the search for extended EWSB sectors 
such as 2HDMs, where multi-lepton final states may arise from many different production and 
decay processes that would evade detection by searches narrowly focused on kinematics or 
resonantly-produced final states of specific topologies. With a potentially sizable multiplicity 
of rare multi-lepton signatures, an extended Higgs sector therefore provides an excellent case 
study for the sort of new physics that could first be discovered in an exclusive multi-channel 
multi-lepton search at the LHC. 

Multi-lepton searches are already sensitive to Standard Model Higgs production [10], as 
well as the production of a SM-like Higgs in rare decay modes of states with large production 
cross sections [11]. This suggests that these studies may be particularly amenable to searching 
for evidence of extended Higgs sectors. Theories with two Higgs doublets enjoy all of the multi- 
lepton final states available to the Standard Model Higgs, albeit with modified cross sections, 
as well as the multi-lepton final states of additional scalars and cascade decays between scalars 
that often feature on-shell W and Z bosons in the final state. These additional particles give 
rise to numerous new production mechanisms for multi-lepton final states. 

The goal of this paper is to perform a detailed survey of the multi-lepton signals that 
arise in some representative 2HDM parameter spaces. In particular, we will consider four 
different CP- and flavor-conserving 2HDM benchmark mass spectra that have qualitatively 
distinct production and decay channels. For each mass spectrum, we will consider each of 
the four discrete types of 2HDM tree-level Yukawa couplings between the Higgs doublets and 
the SM fermions that are guaranteed to be free of tree-level flavor changing neutral currents 
(FCNCs). A study of the sensitivity to the myriad rare production and decay processes over a 
grid of points in the parameter spaces defining these sixteen representative 2HDMs using stan- 
dard simulation techniques, while in principle straightforward, is computationally prohibitive. 
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So instead we employ a factorized mapping procedure to go between model parameters and 
signatures [12]. In this procedure the acceptance times efficiency for each individual produc- 
tion and decay topology is independently determined from monte carlo simulation, assuming 
unit values for all branching ratios in the decay topology. The production cross section and 
branching ratios are then calculated externally as functions of model parameters. The total 
cross section times branching ratio into any given final state at any point in parameter space 
is then given by a sum over the production cross section times acceptance and efficiency for 
each topology times a product of the branching ratios at that parameter space point. For the 
study here, we simulate the acceptance times efficiency in 20 exclusive multi-lepton channels 
for 222 independent production and decay topologies that arise in the four benchmark 2HDM 
spectra. For each benchmark spectrum we combine the 20 exclusive multi-lepton channels to 
obtain an overall sensitivity as a function of two-dimensional mixing angle parameter spaces 
that characterize each of the four discrete types of flavor conserving 2HDMs. With this, we 
identify regions of 2HDM parameter space that are excluded by the existing CMS multi-lepton 
search [9], as well as those regions where future multi-lepton searches at the LHC will have 
sensitivity. 

Beyond requiring CP-conservation and no direct tree-level flavor violation in the Higgs 
sector, we will not address constraints imposed by low energy precision flavor measurements 
on the parameter space of 2HDMs (see [6] and references therein, and [13] for a very re- 
cent analysis). In general, contributions to loop-induced flavor changing processes, such as 
B — > Xs'j, may be reduced by destructive interference among different loops, so that new 
physics outside of our low-energy effective theory can relax flavor bounds on the 2HDM sector. 
Additionally, with the assumptions employed here, flavor constraints are driven by the mass of 
the charged Higgs, which typically does not play a significant role in the production of multi- 
lepton final states. For the benchmark spectra we consider, the charged Higgs may generally 
be decoupled in mass without substantially altering the phenomenology. More generally, we 
emphasize that our benchmark spectra are intended to qualitatively illustrate the relevant 
topologies for producing multi-lepton final states. Various scalar masses may be raised to ac- 
commodate flavor physics without changing the qualitative multi-lepton signatures, though 
of course particular numeric limits will be altered. 

The outline of the paper is as follows: In section 2, we will briefly review the relevant 
aspects of 2HDMs and define the parameter space within which we will conduct our survey. 
In section 3, we will give an overview of the most interesting production and decay channels 
for 2HDM collider phenomenology which result in multi-lepton final states. Additionally, we 
select benchmark spectra that have a representative set of multi-lepton production and decay 
topologies. Section 4 is devoted to summarizing the multi-lepton search strategy and the 
simulation methods we use. The results of our study are displayed in section 5 where we 
identify the regions of parameter space that are excluded on the basis of the existing CMS 
multi-lepton search with 5 fb _1 of 7 TeV proton- proton collisions [9] as well as those regions 
to which future searches will have sensitivity. In section 6 we suggest some refinements to 
future multi-lepton searches that could enhance the sensitivity to extended Higgs sectors. 
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2 Two Higgs Doublet Models 



The physically relevant parameter space specifying the most general 2HDM is large (for a 
review of general 2HDMs see, for example, [5] and [6]). The goal here is not to consider 
the most general theory, but rather to define a manageable parameter space in which to 
characterize multi-lepton signals. The couplings of physical Higgs states that are relevant to 
the production and decay topologies considered below include those of a single Higgs boson 
to two fermions or two gauge bosons, couplings of two Higgs bosons to a single gauge boson, 
and couplings of three Higgs bosons. Other higher multiplicity couplings do not appear in 
the simplest topologies. 

For simplicity we consider CP-conserving 2HDMs that are automatically free of tree-level 
flavor changing neutral currents. With these assumptions, the renormalizable couplings of a 
single physical Higgs boson to pairs of fermions or gauge bosons, and of two Higgs bosons to a 
gauge boson, are completely specified in terms of two mixing angles, as detailed below. With 
a mild restriction to renormalizable potentials of a certain class described below, couplings 
involving three Higgs bosons are specified in terms of Higgs masses and these same mixing 
angles. 

The absence of tree-level flavor changing neutral currents in multi-Higgs theories is guar- 
anteed by the Glashow- Weinberg condition [14] which postulates that all fermions of a given 
gauge representation receive mass through renormalizable Yukawa couplings to a single Higgs 
doublet. With this condition, tree-level couplings of neutral Higgs bosons are diagonal in the 
mass basis. In the case of two Higgs doublets with Yukawa couplings 

- Vy Ukaw a = Yl (QHiVi* + QHiyfd + LH iy fe + h.c.) (2.1) 

i=l,2 

the Glashow- Weinberg condition is satisfied by precisely four discrete types of 2HDMs dis- 
tinguished by the possible assignments of fermion couplings with either yf = or = 
for each of F = u, d, e. Under this restriction, we can always denote the Higgs doublet that 
couples to the up- type quarks as H u . Having fixed this, we have two binary choices for 
whether the down-type quarks and the leptons in (2.1) couple to H u or H^. Of these four 
possibilities, "Type I" is commonly referred to as the fermi-phobic Higgs model in the limit of 
zero mixing, as all fermions couple to one doublet and the scalar modes of the second doublet 
couple to vector bosons only. "Type II" is MSSM-like, since this is the only choice of charge 
assignments consistent with a holomorphic superpotential. "Type III" is often referred to as 
"lepton-specific," since it assigns one Higgs doublet solely to leptons. Finally, "Type IV" is 
also known as "flipped," since the leptons have a "flipped" coupling relative to Type II. These 
possible couplings are illustrated in Table 1. We will restrict ourselves to these four choices 
as they exhaust all possibilities where tree-level FCNCs are automatically forbidden. 

For any of the CP-conserving 2HDMs satisfying the Glashow- Weinberg condition, the 
coefficient of the couplings of a single physical Higgs boson to fermion pairs through the 
Yukawa couplings (2.1) depend on the fermion mass, the ratio of the Higgs expectation values, 
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2HDM I 


2HDM II 


2HDM III 


2HDM IV 


u 


H u 


H u 


H u 


H u 


d 


H u 


H d 


H u 


H d 


e 


H u 


H d 


H d 


H u 



Table 1. The four discrete types of 2HDM H u and Yukawa couplings to right-handed quarks and 
leptons that satisfy the Glashow- Weinberg condition. By convention H u is taken to couple to right 
handed up-type quarks, and the assignments of the remaining couplings are indicated. 



conventionally denned as tan/3 = (H u )/{H d ), and the mixing angle a that diagonalizes the 
2x2 neutral scalar h—H mass squared matrix. The parametric dependences of these couplings 
on a and f3 relative to coupling of the Standard Model Higgs boson with a single Higgs doublet 
are given in Table 2. The parametric dependence of the couplings of the charged scalar, H^, 
are the same as those of the pseudo-scalar, A. 

The renormalizable couplings of a single physical Higgs boson to two gauge bosons are 
fixed by gauge invariance in terms of the mixing angles in any CP-conserving 2HDM as 

g h w = sin(/3 - a)g v g H w = cos(/3 - a)g v g A w = 9h±w^z = (2.2) 

where for V = W, Z the Standard Model Higgs couplings are gw = g and gz = g/ cos 6w, 
where g is the SU(2)l gauge coupling and Ow the weak mixing angle. The renormalizable 
couplings of two physical Higgs bosons to a single gauge boson are likewise fixed in any 
CP-conserving 2HDM as 

9hZA = \gz cos(/3 - a) g H ZA = ~\dz sin(/3 - a 

9hW^H± = T^gcos{/3 - a) g H w*H± = ±^9 ^in(/3 - a) gAW*H± = 7,9 (2-3) 



2* 

i 1 



2 J v ' 2 - 

None of these couplings involve additional assumptions about the form of the full non- 
renormalizable scalar potential, beyond CP conservation. 

The couplings between three physical Higgs bosons depends on details of the Higgs scalar 
potential. Specifying these therefore requires additional assumptions to completely specify 
the branching ratios that appear in some of the decay topologies discussed below. The main 
goal here is to present multi-lepton sensitivities to 2HDMs in relatively simple, manageable 
parameter spaces. A straightforward condition that fulfills this requirement is to consider 
2HDM Higgs potentials that, in additional to being CP-conserving, are renormalizable and 
restricted by a (discrete) Peccei-Quinn symmetry that forbids terms with an odd number of 
H u or H d fields. The most general potential of this type is given by 

Scalar = ^tft^ + m\R\li d + ^-X-y^Jl^f + ^A 2 (i?J^ d ) 2 + A 3 (^^ u )(^^) 



M{HlH d ){H\H u ) + 



l -\ b (HiH d f + h.c. 



(2.4) 
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Table 2. Tree-level couplings of the neutral Higgs bosons to up- and down-type quarks, leptons, and 
massive gauge bosons in the four types of 2HDM models relative to the SM Higgs boson couplings 
as functions of a and /3. The coefficients of the couplings of the charged scalar H ± , are the same as 
those of the pseudo-scalar, A 



This potential has seven free parameters, which may be exchanged for the overall Higgs 
expectation value, the four physical masses m/j, mu^m^ and m#±, and the two mixing 
angles, a. and f3. So all the Higgs boson couplings in a renormalizable 2HDM with the potential 

(2.4) are, for a given mass spectrum, specified entirely in terms of the mixing angles a and 
13. The couplings of three physical Higgs bosons from the potential (2.4) that are relevant to 
the production and decay topologies studied below are 

9Hhh = + 2mi) cos(/3 — a) (sin 2a/ sin 2/3) 

v " ' 

9HAA = ~ ( cos P c °t P sin a + sin f3 tan f3 cos a) + 2m\ cos(/3 — a)) 
v 

Qhh+h- = ~ ( m H ( cos P c °t & sin a + sin (3 tan f3 cos a) + 2m 2 H± cos(/3 — a)) (2.5) 

We emphasize that the choice of the potential (2.4) is illustrative to allow a simple presen- 
tation in terms of a two-dimensional parameter space of mixing angles for a given physical 
spectrum. Although there is additional parametric freedom available in the most general CP- 
conserving 2HDM potential, the phenomenology is qualitatively similar. The only important 
generalization in the production and decay topologies studied below for the most general CP- 
and flavor-conserving 2HDMs as compared with the assumptions outlined here is that the 
partial decay widths of the CP-even heavy Higgs boson, H, to pairs of lighter Higgs bosons 
become free parameters, rather than being specified in terms of a and f3 through the couplings 

(2.5) . 
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3 Multi-lepton Signals of Two Higgs Doublet Models 



The wide range of possibilities for Higgs boson mass spectrum hierarchies and branching 
ratios in 2HDMs yields a diversity of production and decay channels that are relevant for 
multi-lepton signatures at the LHC. Multi-lepton final states become especially important 
when the decay of one Higgs scalar to a pair of Higgs scalars or a Higgs scalar and a vector 
boson is possible. Of course, the availability of these inter-scalar decays comes at a price, as 
the decaying Higgs must be sufficiently heavy for the decay modes to be kinematically open, 
so that the production cross section is reduced. Performing a full multi-dimensional scan of 
the mass spectra of 2HDMs is not only computationally untenable, but also unnecessary for 
our purposes; most of the salient features may be illustrated by exploring a few benchmark 
scenarios in which all the relevant types of cascade decays are realized. We will focus on 
four such mass spectra with various orderings of the scalar mass spectrum, fixing the lightest 
CP-even Higgs mass at 125 GeV in each case. 

The various 2HDM production and decay topologies that give rise to multi-lepton signa- 
tures fall into two broad categories: those resulting from the direct production and decay of 
an individual scalar, and those resulting from cascades involving more than one scalar. The 
first category includes the resonant four-lepton signals of the Standard Model-like Higgs h, 
from gluon fusion and vector boson fusion production followed by h — > ZZ* with Z^ — > ££. 
Other resonant and non-resonant multi-lepton signals arise from quark-anti-quark fusion 
production of Wh, Zh, along with tth associated production with t — > Wb, all followed by 
h — > WW*, ZZ*,tt with leptonic decays of (some of the) W — > tv, ZM — > tl and r — > Ivv. 
These modes were studied in depth in [10] to obtain multi-lepton limits on the Standard 
Model Higgs and simple variations. The same modes of production and decay are available to 
the heavy CP-even scalar, H, albeit with reduced production cross sections due to its larger 
mass and mixing suppression of some of its couplings. While the branching fractions of these 
modes depend on the parameters of the theory, their existence is robust and common to all 
benchmark spectra we consider. In contrast, the sole multi-lepton mode involving direct pro- 
duction of the pseudoscalar, A, without cascade decays through other scalars is ttA associated 
production followed by t — > Wb and A — > tt with leptonic decays of (some of the) W — )• tv 
and r —> Ivv. And there are no multi-lepton signals resulting from direct production of the 
charged Higgs, , without cascade decays through other scalars. 

Scalar cascades add a variety of new multi-lepton processes, including production and 
decay modes that contribute to some of the same final states that arise from a Standard Model 
Higgs boson. Processes of this type include gluon fusion production of A with A —> Zh, ZH 
followed by h,H WW*,ZZ*,tt with (some of the) W -»■ tv, ZM -> tt, and r -> ivv. 
Another example of this type is gluon fusion and vector boson fusion production of H with 
H -»• AA,hh followed by A -»• tt or h -> bb, WW*, ZZ*, tt with (some of the) W -»• Iv, 
ZM -»• 11 and r — > Ivv. With only a single Higgs doublet, direct Standard Model di-Higgs 
production is a very rare process, but resonant heavy Higgs production and decay into these 
final states can be up to two orders of magnitude larger in 2HDMs. 
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Scalar cascade decays of the heavy Higgs scalar, H, can also contribute to entirely new 
multi-lepton final states that do not arise with a single Higgs doublet. These include gluon 
fusion and vector boson fusion production of H with H — > AA, H + H~ , ZA,WH^ with 
A -> 66, Zh, tt, and if* -> tb, tu, Wh with t -> W6 followed by /i -> 66, WW*, ZZ*,tt with 
(some of the) — > Iv, Z^*' —> and r — >• £z^za These processes can give final states with 
up to six W and/or Z bosons. Similar processes in this same category include gluon fusion 
production of A with A — > ZH followed by H — > hh with h — > bb, WW* , ZZ* ,tt with (some 
of the) W — > tv, Z^*> — > U and r — >• fe'ZA These processes can give final states with up to five 
W and/or Z bosons. 

Direct di-Higgs production of non-Standard Model-like Higgs bosons either with or with- 
out scalar cascade decay processes can also give rise to multi-lepton final states that do 
not arise with a single Higgs doublet. These include quark-anti-quark fusion production 
of Ah, AH, H^A followed by H -»• WW* , ZZ* , tt, AA, and H ± -»• tb,TU,Wh,WA with 
t -)• Wb, and A -> 66, rr, all with h,H ^ WW* ,ZZ* ,tt with (some of the) -> 
^(*) an( j T _). ^j/jy, xhe existence of some of these processes is sensitive to mass hier- 

archies in the Higgs spectrum; other production and decay processes of this type can arise 
depending on mass orderings. 

Additional multi-lepton final states not associated with a single Higgs doublet can arise 
from production of non-Standard Model-like Higgs bosons in association with top quarks. 
These include ttH, ttA, and tbH associated production with t — > Wb followed by H — > AA, 
and H± -> Wh, W A, and A -> Zh, bb, tt, all with h, H -> WW*, ZZ* , tt with (some of the) 
W — > iv, Z^*> —> ££ and r — > Ivv. While the production and decay processes listed here and 
above do not completely exhaust all possibilities for contributions to multi-lepton signatures 
in every conceivable 2HDM mass spectrum, they do include the leading topologies for a very 
wide range of mass hierarchies. 

All of the production and decay processes outlined above are represented in one or more of 
the benchmark Higgs mass spectra described below. The value of the scalar masses chosen for 
each benchmark spectrum are shown in Table 3. In the benchmark spectra 1-3, for simplicity 
the pseudoscalar and the charged Higgs are taken to form an isotriplet with degenerate masses. 
In spectrum 4, this simplifying assumption is relaxed, with the pseudoscalar Higgs taken to be 
the lightest scalar. For all four 2HDM spectra, the light, CP-even scalar, h, has no available 
decay modes beyond those of a Standard Model Higgs boson, although the branching fractions 
may significantly differ from the SM values. 

The simplest benchmark spectrum is that with all the heavy non-Standard Model like 
Higgs bosons decoupled. In this case the remaining Standard Model Higgs boson can be 
produced in gluon fusion, vector boson fusion, and in assocation with vector bosons and top 
quarks, and it can decay to h — )■ WW* , ZZ* ,tt. The leading topologies that contribute to 
multi-lepton signatures from these production and decay channels are given in Table 4. These 
topologies are associated to the Standard Model-like Higgs boson in all 2HDMs. The impor- 
tant additional production and decay channels that contribute to multi-lepton signatures 
(beyond those of the Standard Model-like Higgs boson) in each of our four 2HDM benchmark 
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Table 3. Higgs boson masses in the SM Benchmark and our four 2HDM Benchmark Spectra, 
spectra are as follows: 

Benchmark spectrum 1: The heavy neutral Higgs, H, is produced mainly through 
gluon fusion and vector boson fusion, and can decay through the same channels as a heavy 
Standard Model Higgs, plus the new kinematically allowed decay H — > hh. The pseudoscalar, 
A, is produced mainly through gluon fusion and can decay by A — > Zh, ZH. The charged 
Higgs, H^ , does not play an important role in this spectrum. The complete list of topologies 
that contribute to multi-lepton signatures from these production and decay channels, along 
with those from the Standard Model-like Higgs boson, are given in Table 5. 

Benchmark spectrum 2: This spectrum is qualitatively similar to the first, but with 
H — )• hh no longer kinematically allowed. Production of the Heavy Higgs, H, can proceed 
through gluon fusion, vector boson fusion, and in association with vector bosons and top 
quarks, with decays to Standard Model channels. Production of the pseudoscalar, A, through 
gluon fusion production and in association with top quarks with A — > Zh, ZH, tt is much 
greater than in spectrum 1 due to the lower A mass. The charged Higgs, H^, can also be 
produced in association with a top quark, and can decay by H^ — > Wh. The complete 
list of topologies that contribute to multi-lepton signatures from these production and decay 
channels, along with those from the Standard Model-like Higgs boson, are given in Table 6. 

Benchmark spectrum 3: This spectrum is the most rich in the multiplicity of multi- 
lepton final states, as the decay channels H — > hh, AA, H + H~ , AZ are all kinematically open, 
in addition to the Standard Model decay channels. The heavy Higgs, H, can be produced 
in gluon fusion and vector boson fusion. The pseudoscalar, A, is produced in gluon fusion, 
as well as from decays of the H, with decays A — >• Zh,rr. The charged Higgs, H^ , can be 
produced in association with a top quark, or from decay of H with decays H^ — > tv, Wh. 
This spectrum includes topologies with sequential cascade decays through up to three Higgs 
scalars. The complete list of topologies that contribute to multi-lepton signatures from all 
these production and decay channels, along with those from the Standard Model-like Higgs 
boson, are given in Table 7. 

Benchmark spectrum 4: This spectrum breaks the degeneracy between the pseu- 
doscalar, A, and the charged Higgs, H^, in order to highlight the role of a light pseudoscalar. 
Quark-anti-quark fusion production of A with the scalar Higgses, H, h or charged Higgs, 
is significant, with decays A — >■ bb, tt and H^ — > tv, Wh, W A as well as H — > AA, in 
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addition to the Standard Model decay channels. The later decay yields a topology with three 
pseudoscalar Higgses in the final state. The pseudoscalar, A, as well as H and H^, can also 
be produced in association with top quarks. The heavy Higgs, H, can also be produced in 
gluon fusion and vector boson fusion. The very small partial width for the decay h — > AA* in 
this spectrum will be ignored. The complete list of topologies that contribute to multi-lepton 
signatures from all these production and decay channels, along with those from the Standard 
Model-like Higgs boson, are given in Table 8. 

All 233 production and decay topologies listed in Tables 4-8 were individually simulated 
in our studies of multi-lepton signatures of the Standard Model Higgs and our four 2HDM 
spectra benchmarks. Certain channels for the 2HDM benchmarks were omitted for the sake 
of conciseness. In general, channels were omitted if the production cross section times fixed 
Standard Model branching ratios to multi-lepton final states was much less than 1 fb even in 
the most promising regions of parameter space. For nominal simplicity, for the 2HDM bench- 
marks, we omitted associated production channels for h with h — > ZZ* , having found in [10] 
that with the integrated luminosity considered here, these channels did not contribute signif- 
icantly to even low-background search channels. However, with significantly more integrated 
luminosity these channels would begin to contribute to the sensitivity. 



Production 


Decay 


gg-> h 

VBF -> h 
qq -» Wh 
qq — > Zh 
tth 


h^M 

Wh WWW, WZZ, Wtt 
Zh ZWW, ZZZ, Ztt 
tth tiWW, tiZZ, Utt 



Table 4. The 11 independent production and decay topologies simulated for the Standard Model 
Higgs Boson with wih = 125 GeV. The Higgs boson branching ratios are factored out of each topology. 
All top-quark, r-lepton, and W- and Z bosons branching ratios are Standard Model. 



4 Search Strategy and Simulation Tools 

In principle, it might be possible to design a multi-lepton search with sensitivity specifically 
tailored to certain features of the signatures that arise from some of the production and decay 
topologies of 2HDMs. However, designing such a dedicated search would require a detailed 
understanding of backgrounds in many channels that is well beyond the scope of a theory-level 
study. Instead, as done previously in a study of the multi-lepton signatures of the Standard 
Model Higgs boson [10], we will adopt the selection cuts and background estimates of an 
existing CMS multi-lepton analysis [8, 9] to demonstrate the efficacy of a 2HDM multi-lepton 
search. In the conclusions, we will comment briefly on how a focussed search could be further 
optimized to maximize sensitivity to multi-lepton final states arising from an extended scalar 
sector. 
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Production 


Decay 


99 -> h 


h-> 


4^ 




VBF -> h 


h-> 


4^ 




99 


i? - 


>• 4^ 






i? - 


► ft/i - 


> 4VF, WWtt, 4r, ZZ66, ZZ1F W, 4Z, ZZtt 


VBF -> i? 


i? - 


> 4^ 






i? - 


> hh - 


> 4VF, WWtt, At, ZZbb, ZZWW, AZ, ZZtt 




A^ 


Z/i- 


> ZWW, Ztt, ZZZ 




A^ 


Z# - 


-> ZWW, Ztt, ZZZ 




A -> 


ZH - 


-> Zhh -> ZWWWW, ZWWtt, Ztttt, ZZZbb, ZZZWW, 5Z, ZZZtt 


qq -> Wft. 


Wft 


-> WWW, Wtt 


-> Z/l 


Z/i- 


-> ZWW, Ztt 


tth 


tt/i -> ttW W, ttrr 



Table 5. The 37 independent production and decay topologies simulated for the 2HDM Benchmark 
Spectrum 1 with m,h = 125 GeV, ran = 300 GeV, tua = m H ± = 500 GcV. All Higgs boson branching 
ratios are factored out of each topology. All top-quark, 6-quark, r-lepton, and W- and Z-boson 
branching ratios are Standard Model. 



Production 


Decay 


99 -> h 




Al 


VBF -> h 




Al 


99 -> H 


i? -> 


Al 


VBF H 


# -> 


Al 


99 -> -4 


A^ 


Zh -> ZWW, Ztt, ZZZ 




A -> 


Z# ZWW, Ztt, ZZZ 


-> W/i 




-> WWW, Wtt 


-4- Z/i 


Z/i -> ZWW, Ztt 


<7<7 — > WiT 




-> WWW, Wtt 


Zif 


ZH 


-> ZWW, Ztt 


tth 


tth - 


* ttVFPF, ttTT 


ttH 


ttH -> ttWW, ttrr 


ttA 


ttA - 


-> ttTT 




ttA - 


4- ttZ/i ttZVFVF, ttZTT, ttZ&fc, ttZZZ 




ttA - 


-» ttZi? -> ttZVFVF, ttZTT, ttZbb, ttZZZ 


tbH± 


tfcif ± -> tbWh -> tbWWW, tbWTT, tbWZZ 



Table 6. The 34 independent production and decay topologies simulated for the 2HDM Benchmark 
Spectrum 2 with rrih = 125 GeV, ran — 140 GeV, tua = m H ± = 250 GeV. All Higgs boson branching 
ratios are factored out of each topology. All top-quark, 6-quark, T-lepton, and W- and Z-boson 
branching ratios are Standard Model. 

Although the CMS analysis includes hadronically decaying T-leptons, for simplicity of 
simulation, we will consider only strictly leptonic i = e, fi final states (of course, still including 
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Production 


Decay 


99 -> h 


ft- 


-> 4^ 


VBF -> ft 


ft- 


-> 4£ 


99^ H 


i? 


-> 4^ 






-^hh^ AW, WWtt, At, ZZbb, ZZWW, AZ, ZZtt 






-» AA -> 4t 






-> AA -> rrZft -> ttZVFVF, ttZtt, ttZ66, ttZZZ 






AA ^ ZhZh -> ZZWWWW, ZZWWtt, ZZWWbb, ZZttU), ZZtttt 






-> AA -> ZftZft -> ZZ&fote, ZZZZbb, ZZZZtt, ZZZZWW, 6Z 




H 


-> iJ+iJ- -> IWftVFft -> WWWWWW, WWWWtt, WWWWbb, WWtttt 




H 


-> H+R- -> WhWh -> WWrrbb, WW ZZbb, WWWWZZ, WWZZZZ, WW ZZtt 




H 


-> iJ+iJ- -> tz/W/i ™FFVFVF, t^VFtt, tjAFZZ 




H 


-> if -> fW/i -> tbWWW, i&W^TT, iMFZZ 




H 


-> ZA -> Ztt 




H 


-> ZA ZZft -> ZZtt, ZZWW, ZZ66, ZZZZ 




H 


-> Wif ± -> WWft -> VFVFtt, WWWW, WWZZ 


VBF -> i? 


H 


-> 4^ 




H 


-^hh^ AW, WWtt, At, ZZbb, ZZWW, AZ, ZZtt 




H 


-> AA -> 4t 




H 


->• AA ttZ/i -> ttZVFVF, ttZtt, ttZ66, ttZZZ 




H 


-> AA -> ZftZft -> ZZWWWW, ZZWWtt, ZZWWbb, ZZttM>, ZZtttt 




H 


-> AA ZftZft -> ZZbbbb, ZZZZbb, ZZZZtt, ZZZZWW, 6Z 




H 


-> iJ+iJ- -> WftWft -> WWWWWW, WWWWtt, WWWWbb, WWtttt 




H 


-> -> VFftVFft -> WWtt&&, WW ZZbb, WWWWZZ, WWZZZZ, WW ZZtt 




H 


-> H+H- -> tz/W/i ™FFVFVF, t^VFtt, Ti/ffZZ 




H 


-> if -> fW/i -> tbWWW, tfeW^TT, iMFZZ 




H 


-> ZA Ztt 




H 


-> ZA ZZft ZZtt, ZZWW, ZZbb, ZZZZ 




H 


-> VFi?* -> T^T4^ft -> VFVFtt, WWWW, WWZZ 




A 


-> Zft ^ ZT^W^, Ztt, ZZZ 


-> Wft 


Wh -> WWW, Wtt 


q<7 -> Zft 


Zh -> ZVFVF, Ztt 


tth 


ttft -> ttVFW, ttTT 


ttA 


ttA -> ttrr 




ttA — > ttZft ttZVFVF, WZtt, ttZ66, ttZZZ 


tbB ± 


t&ff -> i&VFft -> tbWWW, ife^rr, ibVFZZ 



Table 7. The 111 independent production and decay topologies simulated for the 2HDM Benchmark 
Spectrum 3 with m/i = 125 GeV, mn = 500 GeV, m,A = wih± = 230 GcV. All Higgs boson branching 
ratios are factored out of each topology. All top-quark, 6-quark, T-lepton, and W- and Z-boson 
branching ratios are Standard Model. 



leptonic r decays). Additionally, we treat all hadronic taus as having failed selection criteria, 
thus being identified as jets. Because of this, some events (mainly those involving 4r final 
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Production 


Decay 


99 ->h 


ft^ 4<? 


VBF^h 


ft-> 4^ 


99 -> H 


-> 4^ 




H -> AA -> 4t 


VBF -> 


-> 4<? 




i? -> AA -> 4t 




Wft — > Wtt 




Zh -> ZVKVK, Ztt 


tih 


ttft -> ttWW, ttTT 


ttH 






tiff -> ttAA Utttt, tirrbb 


ttA 


ttA -> ttTT 


tbB ± 


tbB ± -> ffeVK/i -> tbWWW, tfeW^rr, ibVKZZ 




t&H± -> ibVKA i&Wrr 


qq^rH ± A 


H ± A -> WW* -> Wtt66, 




tf± A -> Wftrr -> WWWtt, W^rrrr, VK66rr, ITZZtt 




H^A — > TVTT, tb~TT 




H ± A -> VKAA — > ITtttt, ITtt&6 


qq -> Aft 


Aft -> ttWW, TTTT, ttZZ 


-> Aif 


AH -> ttIT W, tttt, ttZZ 




A£f AAA -> 6r, ttttU> 



Table 8. The 40 independent production and decay topologies simulated for the 2HDM Benchmark 
Spectrum 4 with m h = 125 GeV, m H = 200 GeV, m A = 80 GeV, m H ± = 250 GcV. All Higgs boson 
branching ratios are factored out of each topology. All top-quark, 6-quark, r-lepton, and W- and 
Z-boson branching ratios are Standard Model. 

states) will be categorized differently than in the CMS analysis. For instance, an event with 
three e/fi and one hadronic r that the CMS analysis would have included in a A£ (with It) 
bin, will instead be included in a 3£ bin in our analysis, potentially with higher Ht due 
to the additional energy of the hadronic r-lepton. While this is a deviation from the exact 
procedure of the CMS analysis, it goes in the conservative direction, as the 4£ with It bins 
have significantly smaller backgrounds than the 3£ with Or bins. Thus, if we could implement 
a satisfactory modeling of hadronic r identification in our study, we would expect our bounds 
to become stronger in regions of parameter space where 4r final states are driving the limits. 
For other final states such as H — > hh — > AW, the impact of this effect on our signal is at the 
few percent level or less. 

4.1 Signal channels 

The prompt irreducible Standard Model backgrounds to multi-lepton searches are small and 
arise predominantly through leptonic decays of W and Z bosons. Such backgrounds may 
therefore be reduced by demanding significant hadronic activity and/or missing energy in the 
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events. Hadronic activity can be quantified by the variable Ht, defined as the scalar sum of 
the transverse energies of all jets passing the preselection cuts. The missing transverse energy 
(MET) is the magnitude of the vector sum of the momenta of all particles in the event. 

In order to make use of Hj< and MET, the CMS analysis of [8, 9] divides events with 
H T > 200 (MET > 50) GeV into a high H T (MET) category, and those with H T < 200 
(MET < 50) GeV into a low H T (MET) category. The HIGH H T and HIGH MET require- 
ments (individually or in combination) lead to a significant reduction in Standard Model 
backgrounds. 1 

Another useful observable in reducing backgrounds is the presence of Z candidates, specif- 
ically the existence of an opposite-sign same- flavor (OSSF) lepton pair with an invariant mass 
between 75 — 105 GeV. Events are thus further subdivided, and assigned a No Z channel if 
no such pair exists. It is also useful to characterize events according to whether they may 
contain off-shell 7* / Z* candidates, given by the number of OSSF lepton pairs. Thus, for 
instance, three- lepton events are assigned to the DY0 (no possible Drell-Yan pairs) or DY1 
category (one OSSF pair). The full combination of 3 and 4 lepton events results in 20 pos- 
sible categories of Ht high/low; MET high/low; Z/no Z; and DY0/DY1. The 20 channels 
are presented in Table 10. For each of the 31 and 4i categories, channels are listed from 
top to bottom in approximately descending order of backgrounds, or equivalently ascending 
order of sensitivity, with the last such channel at the bottom dominated by Standard Model 
backgrounds. Events are entered in the table exclusive-hierarchically from the top to the 
bottom. This ensures that each event appears only once in the table, and in the lowest possi- 
ble background channel consistent with its characteristics. Although the backgrounds in the 
individual channels vary over a wide range, all 20 channels are used to compute sensitivity 
limits. 

4.2 Simulation 

For simulating signal processes, we have used MadGraph v4 [15, 16]. In order to simulate 
a general 2HDM in MadGraph, we treat the 2HDM as a simplified model using a modified 
version of the 2HDM4TC model file [17]. Cascade decays were performed in BRIDGE [18]. 
Subsequent showering and hadronization effects were simulated using Pythia [19]. Detector 
effects and object reconstruction was simulated using PGS [20] with the isolation algorithm 
for muons and taus modified to more accurately reflect the procedure used by the CMS 
collaboration. In particular, we introduce a new output variable called trkiso for each muon 
[21]. The variable trkiso is defined to be the sum pr of all tracks, ECAL, and HCAL deposits 
within an annulus of inner radius 0.03 and outer radius 0.3 in AR surrounding a given muon. 
Isolation requires that for each muon, i=trkiso/pr of the muon be less than 0.15. The 
efficiencies of PGS detector effects were normalized by simulating the mSUGRA benchmark 
studied in [8] and comparing the signal in 31 and 4i channels. To match efficiencies with 

1 In the CMS study, a separate binning is also considered using St, a variable defined to be the scalar sum 
of MET, Ht, and leptonic Pt [8]. For simplicity, we will not make use of St here. 
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the CMS study, we applied a lepton ID efficiency correction of 0.87 per lepton to our signal 
events. As discussed earlier, we applied preselection and analysis cuts in accordance with 
those in [8]. 

In order to assess the multi-lepton signatures of the 2HDMs studied here we employ 
a factorized mapping procedure [12] to go between model parameters and signatures. In 
this procedure the acceptance times efficiency is independently determined in each of the 20 
exclusive multi-lepton channels by monte carlo simulation of each individual production and 
decay topology in each of the four 2HDM mass spectra as well as for the individual topologies 
of the Standard Model Higgs boson. The cross section times branching ratio times acceptance 
and efficiency in any of the 20 exclusive channels at any point in parameter space in a given 
mass spectrum is then given by a sum over the production cross section times acceptance and 
efficiency for each topology of that spectrum, times a product of the branching ratios that 
appear in each topology 



where / is a given exclusive final state channel, t labels the topology, and a the branching 
ratios of the decays in the i-th topology. Dependence on the parameter space characterized 
by a and /3 enters only through the production cross sections and decay branching ratios. 
The factorized terms in (4.1) are determined as follows: 

• Acceptance times Efficiency: For each individual production and decay topology 
listed in Tables 4-8, the acceptance times detector efficiency into each of the 20 
exclusive multi-lepton channels listed in Table 10 was simulated with the monte carlo 
tools described above. The acceptance times efficiency of each topology was calculated 
assuming unit branching ratios for all Higgs boson decays but with Standard Model 
values for decays of W and Z bosons, and top quarks and T-leptons. A total of 50,000 
events were simulated for each topology to ensure good statistical coverage of all the 
exclusive multi-lepton channels. 

• Cross Sections: For the case of the Standard Model Higgs boson, the NLO production 
cross sections for gluon fusion, vector boson fusion, and production in association with 
a vector boson or top quarks are taken from the LHC Higgs Cross Section Group [22]. 
For the 2HDM spectra the ratio of LO production partial widths in each production 
channel for h and H relative to a Standard Model Higgs boson of the same mass are 
calculated analytically from the couplings presented in section 2 as functions of the 
mixing parameters a and (3. The NLO Standard Model Higgs production cross sections 
in each production channel are then rescaled by these factors to obtain an estimate for 
the NLO cross sections; for instance the a, (3 dependent cross section for gluon fusion 
production of H is taken to be 




(4.1) 



t 



a 




(4.2) 



- 15 - 



The same procedure of normalizing to Standard Model Higgs boson NLO cross sections 
through the a and j3 dependent ratios of LO production partial widths is used for 
production of A by gluon fusion or in association with top quarks. This is expected 
to be a good approximation since the fractional size of NLO corrections in these cases 
should not be strongly dependent on the parity of the Higgs scalar. For the modes that 
involve production of two Higgs bosons, or of the charged Higgs in association with a 
top quark, the LO cross sections are calculated using Madgraph v4 with a conservative 
i^-factor of K = 1.2 applied. These cross sections are calculated for a single canonical 
value of a and (3 and then rescaled analytically using the couplings in section 2 to obtain 
the cross sections at general values. 

• Higgs Bosons Branching Ratios: For the case of the Standard Model Higgs boson, 
the NLO partial decay widths and branching ratios are taken from the LHC Higgs Cross 
Section Group [22]. For the 2HDM spectra the ratio of LO partial decay widths for h 
relative to a Standard Model Higgs boson of the same mass are calculated analytically as 
functions of the mixing parameters a and (3 using the couplings presented in section 2. 
The NLO Standard Model Higgs boson partial decay widths are then rescaled by these 
factors to obtain estimates for the NLO partial widths; for instance the a, (3 dependent 
partial width for the light scalar h to bb is taken to be 



The same procedure of normalizing to Standard Model Higgs boson NLO partial decay 
widths through the ratio of LO decay widths is used for the H and A decay modes listed 
in Table 9 that are in common with the h decay modes. This estimate is used since, just 
as for a production cross section, the fractional size of NLO corrections to decay widths 
in these cases should not be strongly dependent on the parity of the Higgs scalar. For 
the remainder of the H and A decay modes listed in Table 9 that are kinematically 
open in a given spectrum, as well as the decay modes given in the Table that are 
open, the LO decay widths are calculated analytically [23] as a function of a and /3 
using the couplings in section 2. Except for the charged Higgs decays to quarks, none 
of these decay modes involve strongly interacting particles, so LO widths should be a 
good approximation in this case. The partial widths for all the open decay modes of 
each Higgs scalar in Table 9 are then used to calculate the a and (3 dependent total 
widths and branching ratios in each mass spectrum. 

Using this factorized mapping procedure, each of the 20 exclusive multidepton channels for 
a given benchmark spectrum over the entire a, (3 plane in all four 2HDM types is covered by 
a single set of monte carlo samples for the production and decay topologies. 

In some cases, particularly in Spectrum 3, the total widths of some scalars (particularly 
H) increase drastically in certain regions of parameter space, typically due to enhanced scalar 
couplings. Our simulation and normalization techniques, however, treat all particles in the 




(4.3) 
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Higgs Boson 


Decay Modes 


h 


bb,cc,TT,WW*,ZZ*,gg, 11 ,Z 1 


H 


tt, bb, cc, tt, WW^\ZZ^*\hh, AA, H + H~, ZA, WH ± ,gg, 77, Z 7 


A 


tt, bb, cc, tt, Zh, ZH, gg, 77, Z7 




tb,ts,cs,Tv, WA,Wh,WH 



Table 9. Decay modes of the Higgs boson scalars used in branching ratio calculations. Partial widths 
of the kincmatically open decay modes are calculated in each benchmark spectrum as a function of 
the mixing parameters a and j3 to determine the total width and individual branching ratios. 

narrow width approximation and assume the validity of perturbation theory in the scalar 
couplings. In the regions of parameter space where scalar widths grow large, one expects 
higher-order effects to modify the limits; in this respect the limits we find in high-width 
regions should be viewed as rough estimates subject to potentially large corrections beyond 
the scope of our approach. 

5 Results 

In this section, we present the results of the analysis outlined above using the CMS multi- 
lepton search based on 5 fb _1 of 7 TeV proton- proton collisions at the LHC [9]. We first 
consider the sensitivity of the CMS multi-lepton search to a Standard Model Higgs boson near 
125 GeV before presenting limits in the full 2HDM parameter space for our four benchmark 
spectra. 

For each benchmark, we briefly discuss the major processes that contribute to multi- 
lepton final states, including direct production and decay of individual scalars as well as 
cascades among scalars. We also illustrate many of the partial widths and a ■ Br's for key 
scalar cascades, which helps to capture the qualitative shape of the multi-lepton limits in 
the space of (sina, tan/3). In many cases, the signals of Type I and Type III 2HDM (and 
separately Type II and Type IV 2HDM) are often similar, up to final states involving r- 
leptons. These similarities arise because in each case the quark couplings are identical for 
the pairs of 2HDM types, so in particular the scaling of the h — >■ bb partial widths that often 
govern the total width (as well as the htt couplings that governs the gluon fusion production 
rate) are identical. The only substantial distinction arises in standard channels with r final 
states, since the lepton couplings differ among these pairs of 2HDM types. 

In each case, we show the regions of parameter space excluded by the 5 fb" 1 CMS multi- 
lepton search. In regions not yet excluded, we show the 95% CL limits on the production cross 
section times branching ratio in multiples of the theory cross section times branching ratio for 
the benchmark spectrum and 2HDM type. To compute our 95% CL limits, we used a Bayesian 
likelihood function assuming poisson distributions for each of the 20 channels with a flat 
prior for the signal. We treated the magnitude of the backgrounds in each exclusive channel 
as nuisance parameters with distributions given by a truncated positive definite Gaussian 
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distribution with width equal to the background uncertainty. The number of signal events 
in each exclusive channel for a given a and f3 was obtained from the cross section times 
branching times acceptance and efficiency in each channel times the integrated luminosity. 
For simplicity, we assumed there was no error on the signal. To generate the expected limits, 
a large number of background-only pseudo-experiments were used in place of data. 

For comparison, we also show regions where the heavy, CP-even scalar, H, is currently 
excluded by standard Higgs searches at 7 TeV [1] at roughly the same luminosity of the 
multi-lepton search. For Spectra 1, 3, and 4 we use the combined CMS Higgs limit at 5 fb _1 
of 7 TeV collisions, which is driven by ZZ and WW final states. For Spectrum 2, where 
m H = 140 GeV, we use the WW 2£2u CMS Higgs limit at 5 fb" 1 of 7 TeV collisions, 
which dominates the exclusion limit at this mass. We also consider direct limits on the 
pseudoscalar A and the charged Higgses H^, but these do not impact the parameter space 
explored here. For the pseudoscalar, the best current CMS limits come from MSSM Higgs 
searches for bbA associated production with A — > tt [24]. For a Type II 2HDM, the current 
exclusion is relevant only for tan (3 > 10, and in all other 2HDM types the a ■ Br for bbA 
associated production with A — > tt is smaller than in the Type II case. Searches for di-tau 
resonances [25] do not lead to meaningful limits. Finally, searches for charged Higgses such 
as [26] are sensitive only to production in decays of the top quark, which are not relevant 
for the benchmark spectra considered here. 

5.1 Standard Model Higgs 

We begin by briefly considering the multi-lepton signals of a Standard Model Higgs boson. 
This is useful both as an update to the multi-lepton Higgs search proposed in [10] and as a way 
of understanding certain aspects of the 2HDM multi-lepton signals. In the alignment limit 
defined by sin(/3 — a) = 1 the Higgs expectation values and physical CP-even h eigenstate 
are aligned, and the tree-level couplings of h are identical to those of the Standard Model 
Higgs boson. So in the alignment limit, a 2HDM has an irreducible contribution to multi- 
lepton signatures that is equal to that of the Standard Model Higgs boson, with additional 
contributions coming from the heavier Higgs bosons. The decoupling limit is a special case 
of the alignment limit in which the heavy Higgs scalars are decoupled with large masses. In 
this respect the Standard Model Higgs multi-lepton signals represents a lower bound over a 
sub-space of the 2HDM parameter space, and a limit of the general spectrum space. 

For the Standard Model Higgs, we consider the resonant channels gg —> h — > ZZ* — > 41 
and qq — > h — > ZZ* — > 4£; the non-resonant channels gg —> h — > ZZ* — > 2£2t and qq — > h — > 
ZZ* — > 2£2t; and the associated production channels Zh, Wh, and tth with h — > ZZ* , WW* , 
and tt, all with many possible states yielding multi-lepton signatures. The combined signal 
expectations for a Higgs at 125 GeV in each of the 20 exclusive multi-lepton channels are 
shown in Table 10. As 3£ bins require exactly 3 leptons and 41 bins require > 4 leptons, each 
event appears in the table only once. Although limits may be placed on the signal from any 
individual channel in the multi-lepton search, the greatest sensitivity comes from combining 
all exclusive channels. Combining all multi-lepton channels, we find that the 5 fb _1 multi- 
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Observed 


Expected 


SM Higgs 
Signal 


4 Leptons 














tMET HIGH 


HT 


HIGH 


No Z 





0.018 ± 0.005 


0.03 


tMET HIGH 


HT 


HIGH 


Z 





0.22 ± 0.05 


0.01 


tAyTTTT 1 1 1 T ^ 1 1 
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1 1 1 
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U.Uu 


tMET HIGH 


HT 


LOW 
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UT 
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1 1 If M 1 

1 1 1 Vj 1 1 




n 
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U.UUD it U.UUl 
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tMET LOW 


HT 


HIGH 


z 


i 


0.83 ± 0.33 


0.01 


tMET LOW 


HT 


LOW 


No Z 


i 


2.6 ± 1.1 


0.36 


tMET LOW 


HT 


LOW 


Z 


33 


37 ± 15 


1.2 


3 Leptons 














tMET HIGH 


HT 


HIGH 


DYO 


2 


1.5 ± 0.5 


0.15 


tMET HIGH 


HT 


LOW 


DYO 


7 


6.6 ± 2.3 


0.67 


f MET LOW 


HT 


HIGH 


DYO 


1 


1.2 ± 0.7 


0.04 


tMET LOW 


HT 


LOW 


DYO 


14 


11.7 ± 3.6 


0.63 


tMET HIGH 


HT 


HIGH 


DY1 No Z 


8 


5.0 ± 1.3 


0.38 


tMET HIGH 


HT 


HIGH 


DY1 Z 


20 


18.9 ± 6.4 


0.19 


tMET HIGH 


HT 


LOW 


DY1 No Z 


30 


27.0 ± 7.6 


1.8 


MET HIGH 


HT 


LOW 


DY1 Z 


141 


134 ± 50 


1.6 


tMET LOW 


HT 


HIGH 


DY1 No Z 


11 


4.5 ± 1.5 


0.13 


tMET LOW 


HT 


HIGH 


DY1 Z 


15 


19.2 ± 4.8 


0.09 


MET LOW 


HT 


LOW 


DY1 No Z 


123 


144 ± 36 


1.8 


MET LOW 


HT 


LOW 


DY1 Z 


657 


764 ± 183 


4.3 



Table 10. Observed and expected number of events in various exclusive multi-lepton channels from 
the CMS multi-lepton search with 5 fb" 1 of 7 TeV proton-proton collisions [9], along with expected 
number of Standard Model Higgs boson signal events for irih — 125 GeV after acceptance and efficiency. 
HIGH and LOW for MET and HT indicate $ T > 50 GeV and H T > 200 GeV respectively. DYO 
= f±«T ; DY1 = e ± e+e-,e ,± £+£-, for i = e,^. No Z and Z indicate \m u - m z \ > 15 GeV for any 
opposite sign same flavor pair. The channels with moderate to good sensitivity to multi-lepton Higgs 
boson signals arc indicated with daggers. 

lepton CMS results [9] yield the expected and observed limits for a Standard Model Higgs at 
nih = 120, 125, and 130 GeV shown in Table 11. The dominant decay modes and exclusive 
channels contributing to these limits were discussed in detail in [10]. 

The multi-lepton signals of h remain important in the general 2HDM parameter space, 
both through Standard Model production of h and the production of h in scalar cascades. 
The variation in these signals as a function of sin a and tan f3 for the four types of 2HDM was 
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120 GeV 


125 GeV 


130 GeV 


Observed 


5.4 


4.9 


3.5 


Expected 


4.2 


3.8 


2.8 



Table 11. Observed and expected 95% CL limits from the CMS multi-lepton search with 5 fb _1 of 
7 TeV proton-proton collisions [9] on the Higgs boson production cross section times branching ratio 
in multiples of that for Standard Model Higgs multi-lepton production and decay topologies listed in 
Table 4 with Standard Model branching ratios. Limits are obtained from an exclusive combination of 
the observed and expected number of events in all the multi-lepton channels presented in Table 10. 

studied in detail in [27]; in what follows, we will often refer to these results to understand the 
parametric changes in the multi-lepton limit across the 2HDM parameter space. 

5.2 Spectrum 1 

Now let us turn to the multidepton signals and limits of our 2HDM benchmark spectra. The 
multidepton limits on the first benchmark spectrum for all four types of 2HDM are shown in 
Figure 1. Limits in this and the following figures were obtained from an exclusive combination 
of the observed and expected number of events in all the multidepton channels presented in 
Table 10 on an evenly-spaced grid in — 1 < sin a < and 1 < tan/3 < 10 with spacing 
A(sina) = 0.1 and A(tan/3) = 1; contours were determined by numerical interpolation 
between these points. 

In addition to the Standard Modeldike production and decays of scalars to SM final 
states, the first benchmark spectrum also features the inter-scalar decays H — > hh, A —> Zh, 
and A —> ZH. The partial widths for these three inter-scalar decays (which are independent 
of the 2HDM type) and the a ■ Br for the dominant processes gg —> H — > hh, gg —> A — )■ Zh, 
and gg — > A — > ZH (which depend weakly on the 2HDM type; here, we display those of a 
Type I 2HDM) are shown in Figure 2; their parametric behavior as a function of sin a and 
tan/3 helps to explain many of the detailed features of the exclusion limits in Figure 1. 

The partial width, T(H — > hh), has a complicated dependence on a, [3, but is greatest 
when tan/3 is large and sin a ~ —0.85. This process only contributes significantly to multi- 
lepton limits in 2HDM types for which the multi-lepton decays of h are unsuppressed in the 
same region where Br(if — > hh) is large. The partial width, F(A — > Zh) oc cos 2 (/3 — a), is 
largest away from the alignment limit, while the partial width, T(A — > ZH) oc sin 2 (/3 — a), 
is largest in the alignment limit. In both cases, the multi-lepton limits are strongest for 
2HDM types where the multi-lepton decays of h and H are significant when Br(j4 — > Zh) 
and Br(A — > ZH) are respectively large. 

On the production side, the dominant production cross section for H, a(gg — > H), is 
largest at small tan/3 and sin a — > — 1, while the dominant cross section for A, a(gg — > A), is 
independent of sin a (since the pseudoscalar couplings to fermions, and hence gluons, depend 
only on tan /3) and increases as tan /3 — > 0. These production cross sections and scalar partial 
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Figure 1. Multi-lepton limits from the CMS multi-lepton search with 5 fb _1 of 7 TeV proton-proton 
collisions [9] for the production and decay topologies of Benchmark Spectrum 1 given in Table 5, for 
Type I (top left), Type II (top right), Type III (bottom left), and Type IV (bottom right) couplings 
as a function of sin a and tan/3. Limits were obtained from an exclusive combination of the observed 
and expected number of events in all the multi-lepton channels presented in Table 10. The solid and 
dashed lines correspond to the observed and expected 95% CL limits on the production cross section 
times branching ratio in multiples of the theory cross section times branching ratio for the benchmark 
spectrum and 2HDM type. The blue shaded regions denote excluded parameter space. The solid 
red line denotes the alignment limit sin(/3 — a) = 1. The gray shaded region corresponds to areas of 
parameter space where vector decays of the heavy CP-even Higgs, H — > VV, are excluded at 95% CL 
by the SM Higgs searches at 7 TeV [1]. 

widths are largely independent of the 2HDM type; the gluon fusion rates for Type II and 
Type IV 2HDM increase slightly at large tan j3 due to the sizable bottom quark coupling. 
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Figure 2. 2HDM Benchmark Spectrum 1 partial widths T(H -> hh), T(A -> Zh), and T(A -> ZH) 
in units of GeV, and cross section times branching ratios a ■ Br(gg — >• H — > hh), a ■ Br(gg — >• A — > Zh), 
and a ■ Br(gg — >• A — > Zi/) in units of pb, all for Type I couplings. These partial widths and 
a ■ Brs are qualitatively similar for the other types of 2HDM couplings; the production cross sections 
a(gg —> H,A) are moderately enhanced at large tan/3 for Type II and Type IV 2HDM due to the 
contribution from bottom loops. 



The threefold combination of production rates, inter-scalar decay widths, and multi- 
lepton widths of scalars determines the shape of limits in the plane of sin a and tan/3. These 
vary among different 2HDM types, though similarities between Type I & III and between 
Type II & IV make it worthwhile to discuss these two sets together. 

Types I & III 

In the Type I 2HDM, the multi-lepton signals of the SM-like Higgs, h, generally decrease 
as we move away from the alignment limit (in large part because the coupling to vectors 
is suppressed, reducing both the Vh associated production rate and the branching ratios, 
Br(/i — > VV*); for an extended discussion, see [27]), but are not a strong function of sin a 
and tan/3; only near sin a — > — 1 are the a ■ Br for the conventional multi-lepton channels of 
h significantly diminished. However, the SM-like multi-lepton signals of h are typically never 
enhanced as we move away from the alignment limit (the exception being a mild enhancement 
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of VBF and Vh associated production with h —> VV* at small tan/3 and sin a — > — 1; see 
[27] for more detail). In the region where the multi-lepton signals of h are diminished, the 
conventional multi-lepton signals of H are correspondingly enhanced since the HVV coupling 
is complementary to the HVV coupling. While for = 300 GeV, the production cross 
section for H is somewhat smaller than that of h, it nonetheless contributes significantly to 
multi-lepton limits near sin a — > — 1 through primarily SM-like production and decay modes. 
Note that the direct decays of the pseudoscalar A never result in more than two leptons, so 
the pseudoscalar contributes to the multi-lepton signal only through scalar cascades and ttA 
associated production. 

In addition to the conventional SM-like production and decay modes of h and H, we must 
also consider the various production channels involving inter-scalar decays. The a ■ T$r(gg —> 
H —> hh) is largest at large tan/3 and sin a ~ —0.8 where guhh is largest. The parametric 
behavior of this a ■ Br, along with the fact that the multi-lepton final states of h in a Type I 
2HDM are only mildly suppressed when a ■ Br(<?g — > H — > hh) is significant, largely explains 
the strengthening of the multi-lepton limit around sin a ~ —0.85. 

For the pseudoscalar, a ■ Br(gg — > A — > Zh) is large away from the alignment limit, but 
decreases at large tan (3 due to the falling gluon fusion rate for A. Similarly, a ■ Br(gg — > 
A — > ZH) is large only at low tan/3, since the branching ratio for A — > ZH is large along 
the alignment line but the gluon fusion rate for A again decreases at large tan /3. Thus, both 
a ■ T5i(gg —tA—t Zh) and a ■ Br(gg — )■ A — >■ ZH) contribute to limit-setting at small tan/3, 
essentially independent of sin a, while a ■ Br(gg — > A — > Zh) also contributes at larger tan /3 
for sin a < —0.5. 

All three scalar decays contribute to setting the strongest limits at small tan /3 (relatively 
insensitive to sin a), while a ■ TSr(gg — > H — > hh) predominantly explains the limits at large 
tan/3 around sin a ~ —0.85. The additional contributions from scalar cascades are exemplified 
in Figure 3, which illustrates the Ht and MET distributions for the sum of multi-lepton events 
at the point (sina = — 0.9, tan (3 = 1.0), distinguished by the initial scalar produced in each 
multi-lepton event. 

The multi-lepton signals in the Type III, or "lepton-specific," model are similar to those of 
the Type I model, since the couplings of the Higgs scalars to quarks and vectors are identical 
for these 2HDM types. The exception is a significant improvement in the limits around 
—0.9 < sina < —0.6 relative to the Type I 2HDM. Here, the branching ratio, Br(/i — > rr), 
is substantially increased over the SM rate and contributes both through SM-like associated 
production of h and production of H — > hh with one or both h decaying to rr. Indeed, 
processes such as Vh associated production with h — > tt are as much as ten times larger 
than the SM rate, with a ■ T$r(Wh — > Wtt) as large as several hundred fb. Scalar cascades 
involving rs are even more important, with a ■ Br(<7<? — > H — > hh — > 4r) as large as several 
pb. The enhancement of T(h — > tt) renders this the 2HDM type most amenable to detection 
by the multi-lepton search, and, in fact, a large region of parameter space is already excluded 
by the CMS multi-lepton search with 5 fb . While some of this region is already excluded 
by conventional searches for h — > tt, there exist regions not constrained by current searches 
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Figure 3. The 2HDM signal transverse hadronic energy distribution (left) and missing transverse 
energy distribution (right) after acceptance and efficiency for 7 TeV proton-proton collisions arising 
from the production and decay topologies of Benchmark Spectrum 1 given in Table 5 with mu = 125 
GeV, m H = 300 GeV, m H ± = m A = 500 GeV, for Type I 2HDM couplings with sin a = -0.9 and 
tan/3 = 1.0. Signal events correspond to those falling in the exclusive three- or four-lepton channels 
labelled with a dagger in Table 10 that have moderate to good sensitivity. The colors indicate the 
initial type of Higgs boson produced. For each color, the lighter shade corresponds to three-lepton 
channels, while the darker shade corresponds to four-lepton channels. The bin size is 40 GeV for Ht 
and 10 GeV for $ T , and in both cases the highest bin includes overflow. 

where the dominant multi-lepton limit comes from scalar cascades. 
Types II & IV 

A very important difference in the phenomenology of the Type II & IV 2HDM compared to 
the preceding description of the Type I & III phenomenology is that the down-type quarks 
now couple to rather than H u , thus the partial width of h — > bb has an entirely different 
parametric dependence. Since this decay mode dominates in the SM-like alignment limit, its 
variation sharply affects the Br's of all other decay modes as well. For instance, the multi- 
lepton signals of the SM-like Higgs h change rapidly as we move away from the alignment 
limit, decreasing sharply with increasing tan /3 above the sin(/3 — a) = 1 line due to the rapidly 
increasing partial width, T(h — > bb), and rising rapidly below sin(/3 — a) = 1 as T(h —> bb) 
drops. Thus at large tan /3 above the alignment line, the multi-lepton signals of h diminish 
rapidly, weakening the limit both from SM-like production of h and from new associated 
production, such as H — > hh. The only exception are multi-lepton signals involving h — > tt, 
since T(h —> Tr)/Y(h —> bb) is fixed in a Type II 2HDM. On the other hand, below the 
alignment line there is an overall enhancement of multi-lepton decays involving h — > VV* 
since the partial width T(h —> bb) drops, leading to an increase in the purely SM-like multi- 
lepton production and decay modes of h. As sin a — > — 1, the direct multi-lepton decays of H 
somewhat compensate for the loss of h signals, but there is a wide region of large tan /3 and 
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moderate sin a where neither h nor H decays significantly to multi-lepton final states; this is 
clearly displayed by the weak limits in the range —0.9 < sin a < —0.2. 

Scalar cascade decays do not significantly help to constrain a Type II 2HDM. While the 
a ■ Br(<7<7 — > H — > hh) is parametrically similar to the Type I 2HDM, in a Type II 2HDM the 
SM-like Higgs h decays predominantly to bb in this region, so this channel does not contribute 
substantially to multi-lepton limits (except for the rare hh — > 4r). Likewise, the contributions 
from a ■ Br(gg — > A — > Zh) at large tan/3 lead to multi-lepton signals only through h — > tt. 

At low tan (3, the direct multi-lepton decays of h are still significant, as are the added 
contributions from H — > hh, A — > Zh, and A —> ZH. The multi-lepton limits on the first 
benchmark spectrum for a Type II 2HDM are strongest at low tan/3, where h decays and 
inter-scalar decays to multi-lepton final states are enhanced; limits at sin a — > —1 come 
predominantly from direct decays of H, while those at sin a — > come from direct decays 
of h. The contributions of the pseudoscalar in this limit are exemplified by Figure 4, which 
illustrates the and MET distributions for the sum of multi-lepton events at the point 
(sin a = — 0.3, tan /3 = 1.0), for which there is a large contribution from A — > Zh, ZH. 




Figure 4. The 2HDM signal transverse hadronic energy distribution (left) and missing transverse 
energy distribution (right) after acceptance and efficiency for 7 TeV proton-proton collisions arising 
from the production and decay topologies of Benchmark Spectrum 1 given in Table 5 with m/, = 125 
GeV, m H = 300 GeV, m H ± = m A = 500 GeV, for Type II 2HDM couplings with sin a = -0.3 and 
tan/3 = 1.0. Signal events correspond to those falling in the exclusive three- or four-lepton channels 
labelled with a dagger in Table 10 that have moderate to good sensitivity. The colors indicate the 
initial type of Higgs boson produced. For each color, the lighter shade corresponds to three-lepton 
channels, while the darker shade corresponds to four-lepton channels. The bin size is 40 GeV for Ht 
and 10 GeV for ^? T , and in both cases the highest bin includes overflow. 

The multi-lepton signals in the Type IV, or "flipped," model are similar to that of the 
Type II model, since the couplings of the Higgs scalars to quarks and vectors are identical for 
these 2HDM types. The notable exception are the reduced limits in the region of moderate 
sin a and large tan (3. This reduction in sensitivity is due to the fact that in a Type IV 2HDM 
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the partial width, T(h — > tt), no longer scales with T(h — > bb), and so in the region where 
r(/i — > bb) is particularly large there are no longer meaningful contributions to multi-lepton 
limits from h — > tt with leptonically decaying rs. In particular, this removes possible multi- 
lepton signals from associated production of h in this region, both through SM associated 
production and scalar cascades. 

5.3 Spectrum 2 

The multi-lepton limits on the second benchmark spectrum are shown in Figure 5. Much like 
the first benchmark spectrum, this spectrum includes the scalar decays A —> Zh and A — > ZH, 
albeit with greater cross sections since txia = 250 GeV in this spectrum. However, the decay 
H — > hh is now kinematically forbidden. Since the parametric behavior of the relevant partial 
widths and a ■ Br's is the same as in the first benchmark up to overall rescalings, we do not 
show them explicitly, but emphasize that the cross sections for production of A and H are 
substantially larger compared to the first benchmark since both A and H are lighter in this 
case. 

Types I & III 

The multi-lepton limits for Type I 2HDM are similar to those of the Type I model for Spectrum 
1, albeit without the contributions from H — > hh. Particularly, the stronger limits around 
sin a ~ —0.85 in Spectrum 1 are absent here, but otherwise the parametric contributions are 
similar. The limits for this spectrum are stronger at small tan/3 because the now lighter A 
has a larger production cross section, a(gg —> A), than in Spectrum 1. Similarly, the limits 
are stronger as sin a — > — 1 since here the direct production and multi-lepton decays of H 
dominate the limit, and the production cross section for H is effectively SM-like in this region 
since mjj = 140 GeV. 

Likewise, the multi-lepton limits for Type III 2HDM are similar to those of the Type III 
model for Spectrum 1, although they again lack the contributions from H — > hh, meaning 
that there is no significant 4r contribution with this spectrum. 

Types II & IV 

Unsurprisingly, the limits for Type II & Type IV 2HDM are similar to the analogous limits 
in Spectrum 1, although somewhat stronger due to the enhanced production cross sections 
for A and H. Note that there is no significant weakening of the limit at large tan f3 and 
moderate sin a compared to Spectrum 1, despite the disappearance of the decay H — > hh. 
This exemplifies the fact that in Type II and Type IV 2HDM, the multi-lepton decays of h 
are suppressed in this range, so the presence or absence of H — > hh does not significantly 
alter the limit. 

5.4 Spectrum 3 

The multi-lepton limits on the third benchmark spectrum for all four types of 2HDM are 
shown in Figure 6. The third benchmark spectrum enjoys a plethora of inter-scalar cascade 
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Figure 5. Multi-lepton limits from the CMS multi-lcpton search with 5 fb _1 of 7 TeV proton-proton 
collisions [9] for the production and decay topologies of Benchmark Spectrum 2 given in Table 5, for 
Type I (top left), Type II (top right), Type III (bottom left), and Type IV (bottom right) couplings 
as a function of sin a and tan/3. Limits were obtained from an exclusive combination of the observed 
and expected number of events in all the multi-lepton channels presented in Table 10. The solid and 
dashed lines correspond to the observed and expected 95% CL limits on the production cross section 
times branching ratio in multiples of the theory cross section times branching ratio for the benchmark 
spectrum and 2HDM type. The blue shaded regions denote excluded parameter space. The solid 
red line denotes the alignment limit sin(/3 — a) = 1. The gray shaded region corresponds to areas of 
parameter space where vector decays of the heavy CP-even Higgs, H — > WW* , are excluded at 95% 
CL by the SM Higgs searches at 7 TeV [1]. 



decays. In particular, the important inter-scalar decays include H — > hh, H — > AA, H — > 
H+H-, H -> ZA, F± ->• W ± h, and A -)• Zh. The fact that H -)• H + H~ , AA, ZA and both 
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— > W^h and A — > Zh are open allows for the possibility of multi-step cascades involving 
three Higgs scalars. Also note that the range of possible decays of H means that the overlap 
of large T(H — > hh) with multi-lepton decays of h is not as important to limit-setting as it 
was in Spectrum 1, since, e.g., H — > AA, ZA with A — > tt may be important even when the 
multi-lepton decays of h are small. However, since H is relatively heavy in this benchmark 
{rriH = 500 GeV), the direct multi-lepton decays of H are less important to limit-setting 
relative to other benchmarks due to the lower production cross section. The partial widths 
and a ■ Br for those processes unique to Spectrum 3 are shown in Figure 7 (the parametric 
dependence of H —> hh and A — > Zh were already shown in Figure 2 and the dependence of 

— > W^h will be shown in Figure 11 when we discuss Spectrum 4). 

The partial widths Y{H — > A A) and T(H —> H + H~) are complicated functions of a and 
/3, but grow as tan/3 increases and sin a goes to zero. The partial widths, Y{H — > ZA) and 
T(H — > H^W^), scale simply as sin 2 (/3 — a), and so is largest in the alignment limit, while 
the partial widths, T(A -)• hZ) and T{H ± -)■ W ± h), scale as cos 2 (/3 — a) and is largest away 
from the alignment limit. 

Note in Figure 7 the partial widths, T(H — > A A) and T(H — > H + H~), grow quite large 
with increasing tan/3, such that the total width of H exceeds its mass for tan/3 > 5 and 
sin a > —0.8. In this regime, both the perturbative expansion in scalar couplings and the 
narrow width approximation break down, and the precise exclusion limit should be treated 
with caution. 

On the production end, as noted earlier the dominant production mode for H, a(gg — > H), 
is largest at small tan f3 and sin a — > — 1. The combination of this dependence and the partial 
widths implies that a ■ Br(gg — > H — > AA) and a ■ Br(gg — > H — > H + H~) are largest 
at moderate sin a, peaking around sin a ~ —0.8 and increasing mildly with tan/3; both 
contribute over a somewhat wider range than gg — > H — )■ hh. In contrast, a ■ Br(g<? — > H — > 
ZA) is largest at low tan/3 and sin a — > —1. 

Types I & III 

The signals of the Type I 2HDM for the third benchmark spectrum are similar to those of 
the first benchmark spectrum, to the extent that they are largely governed by the multi- 
lepton final states of h combined with the scalar decays of H and A. However, in contrast to 
Spectrum 1, here the direct multi-lepton decays of H are less significant in limit-setting since 
the production cross section for run = 500 GeV is considerably smaller. Thus, the limits 
at large tan/3 and sin a — > — 1 coming from direct multi-lepton decays of H are noticeably 
weaker in this case. On the other hand, scalar decays of H contribute meaningfully over a 
wide range in sin a since a ■ Br(gg — )■ H — >■ AA) and a ■ Br(gg — > H — > H + H~ ) change slowly 
as a function of sin a compared to a ■ Br(gg — > H — > hh). 

In the case of processes involving H — > AA, the multi-lepton limits are dominated by the 
decays A — > Zh rather than A — > tt. This is because in a Type I model the Att coupling 
decreases with increasing tan/3, so that the branching ratio Br(A — > tt) is not large in the 
same region as a ■ Br(gg — > H — > AA). In contrast, the branching ratio Br (^4 — > Zh) is large 
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Figure 6. Multi-lepton limits from the CMS multi-lcpton search with 5 fb _1 of 7 TeV proton-proton 
collisions [9] for the production and decay topologies of Benchmark Spectrum 3 given in Table 5, for 
Type I (top left), Type II (top right), Type III (bottom left), and Type IV (bottom right) couplings 
as a function of sin a and tan/3. Limits were obtained from an exclusive combination of the observed 
and expected number of events in all the multi-lepton channels presented in Table 10. The solid and 
dashed lines correspond to the observed and expected 95% CL limits on the production cross section 
times branching ratio in multiples of the theory cross section times branching ratio for the benchmark 
spectrum and 2HDM type. The blue shaded regions denote excluded parameter space. The solid 
red line denotes the alignment limit sin(/3 — a) = 1. The gray shaded region corresponds to areas of 
parameter space where vector decays of the heavy CP-even Higgs, H — > VV* , are excluded at 95% 
CL by the SM Higgs searches at 7 TeV [1]. In all cases, for tan j3 > 5 and sin a > —0.8 the total width 
of H grows comparable to its mass and the precise exclusion limit in this region is subject to large 
theoretical uncertainties, these regions are highlighted in light red. 
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Figure 7. 2HDM Benchmark Spectrum 3 partial widths T(H ->■ AA), T(H H+H~), and T(H ->■ 
ZA) in units of GeV, and cross section times branching ratios a -Br(gg — » H — > AA), a-Bi(gg — > H — > 
H + H~), and a ■ Br(gg —> H — > ZA) in units of pb, all for Type I couplings. These partial widths and 
a ■ Brs are qualitatively similar for the other types of 2HDM; the production cross section a(gg — » H) 
is moderately enhanced at large tan/3 for Type II and Type IV 2HDM due to the contribution from 
bottom loops. 

precisely when J5t(H — > AA) is large, hence H — > AA — > ZhZh contributes substantially to 
the limit at large tan/3 and —0.9 < sin a < —0.4, with a ■ Br(gg — > H — > AA —> ZhZh) 
growing as large as ~ 120 fb in the region of study. 

For processes involving H — > H + H~, the multi-lepton limits always require at least 
one charged Higgs to decay via — > W h, since the other decay modes such as e.g. 
H + — > tb, t + v give at most one lepton. In a Type I model, Br(if =t — > W^h) is sizable when 
Bi(H —7- H + H~) is large, so H — > H + H~ — > W + hW~h is important at large tan/3 in the 
range —0.9 < sin a < —0.5. Processes involving H — > H + H~ with one decay to tb and tv 
are also important at moderate tan (3. 

As in previous cases, gg —> A —> Zh is important at small tan /?, as is gg — > H — > ZA with 
both A — ^ tt and A — > Zh. Various exemplary features of the third benchmark spectrum 
with Type I 2HDM couplings are shown in Figure 8, which illustrates the and MET 
distributions for the sum of multi-lepton events at the point (sin a = — 0.9, tan (3 = 1.0), 
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distinguished by the initial scalar produced in each multi-lepton event. 
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Figure 8. The 2HDM signal transverse hadronic energy distribution (left) and missing transverse 
energy distribution (right) after acceptance and efficiency for 7 TeV proton-proton collisions arising 
from the production and decay topologies of Benchmark Spectrum 3 given in Table 7 with mu = 125 
GeV, m H = 500 GeV, m H ± = m A = 230 GeV, for Type I 2HDM couplings with sin a = -0.9 and 
tan/3 = 1.0. Signal events correspond to those falling in the exclusive three- or four-lepton channels 
labelled with a dagger in Table 10 that have moderate to good sensitivity. The colors indicate the 
initial type of Higgs boson produced. For each color, the lighter shade corresponds to three-lepton 
channels, while the darker shade corresponds to four-lepton channels. The bin size is 40 GeV for Ht 
and 10 GeV for |5 T , and in both cases the highest bin includes overflow. 



The Type III 2HDM shares many of the qualitative features of the Type I 2HDM, albeit 
with additional contributions to multi-lepton signals coming from the fact that the partial 
widths T(/i — > tt) and T(A — > tt) grow with tan/3. So, in addition to the significant signals 
discussed earlier, both H — > hh — > 4r and H — > AA — > 4r are important in the Type III 
2HDM, particularly at moderate sin a and large tan/3 where J5r(H — > hh,AA) are large and 
so too are Br(/i, A — > tt). Taken together, these contributions are still not as great as in 
Spectrum 1 due to the reduced production cross section for H, but nonetheless lead to large 
regions already excluded using the 5 fb _1 data. 

Types II & IV 

As in previous cases, the multi-lepton final states of h decrease rapidly above the alignment 
limit, with the sole exception of h — > tt. Here, the reduced contribution from direct multi- 
lepton decays of H is particularly noticeable, with a substantial weakening of the limit as 
sina — > —1. 

Much as in Spectrum 1 Type II, processes involving H — > hh contribute little to the 
limit, since h has suppressed multi-lepton final states when Br(H — > hh) is large. The decay, 
H — > AA, is somewhat more important, but, as with the Type III model, the contribution to 
multi-leptons comes primarily from A — > tt as opposed to A — > Zh, especially at large tan/3. 
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The Att coupling grows with tan/3 in a Type II 2HDM, but, as before, A — > bb, with the 
same parametric scaling, still dominates the total width of A. Similarly, decays primarily 
to tb and tv at large tan j3, so — > W^h is suppressed in this range and processes involving 
H — > H + H~ do not contribute much to the multi-lepton limits. 

The processes gg — >■ A — >■ Zh and gg —> H —> Z(A — > Zh) are important at small 
tan /3; here the multi-lepton decays of h are enhanced below the alignment line, so that these 
processes contribute significantly to the limit through the direct multi-lepton decays of h. The 
contributions of the pseudoscalar are exemplified by Figure 9, which illustrates the Ht and 
MET distributions for the sum of multi-lepton events at the point (sin a = — 0.2, tan /3 = 1.0), 
for which there is a large contribution from A — > Zh. 
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Figure 9. The 2HDM signal transverse hadronic energy distribution (left) and missing transverse 
energy distribution (right) after acceptance and efficiency for 7 TeV proton-proton collisions arising 
from the production and decay topologies of Benchmark Spectrum 3 given in Table 7 with m/, = 125 
GeV, m H = 500 GeV, m H ± = m A = 230 GeV, for Type II 2HDM couplings with sin a = -0.2 and 
tan/3 = 1.0. Signal events correspond to those falling in the exclusive three- or four-lepton channels 
labelled with a dagger in Table 10 that have moderate to good sensitivity. The colors indicate the 
initial type of Higgs boson produced. For each color, the lighter shade corresponds to three-lepton 
channels, while the darker shade corresponds to four-lepton channels. The bin size is 40 GeV for Ht 
and 10 GeV for $ T , and in both cases the highest bin includes overflow. 



The Type IV 2HDM recapitulates many of the features of the Type II 2HDM, albeit 
without significant contributions from h —> tt or A — > rr at large tan/3. This eliminates 
contributions from, e.g., H — > hh — > 4r and H — > AA — > 4r, so that the multi-lepton 
limits are particularly weak at moderate sin a and large tan /3. As before, the multi-lepton 
decays of h are important below the alignment line, and accumulate extra contributions from 
gg ->■ A — ^ Zh and gg ->■ H -> Z(A -> Zh) at low tan/3. 

5.5 Spectrum 4 

The multi-lepton limits on the first benchmark spectrum for all four types of 2HDM are shown 
in Figure 10. The fourth benchmark spectrum highlights the signals of a light pseudoscalar, 
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both through decays of other scalars and through direct production in association with those 
scalars. Kinematically available inter-scalar decays include H — > AA, — > W h, and 

— > W^A, while interesting associated production processes unique to this benchmark 
include qq — > H^A, qq — > Ah, and qq — > AH through off-shell W and Z bosons. The partial 
widths and a ■ Brs for several of these processes are shown in Figure 11. 

The partial width r(if =t — > W^h) scales as cos 2 (/3 — a) and hence grows away from 
the alignment limit. In contrast, r(// =t — > W^A) is entirely independent of the angles a, (3. 
On the production side, cr(qq — > Ah) oc cos 2 (/3 — a) grows away from the alignment limit, 
while (j(qq — > AH) oc sin 2 (/3 — a) grows as we approach the alignment limit. The production 
cross section a(qq — > H^A) is likewise independent of a, /3 since it scales as the square of the 
H^W^A coupling. However, the partial widths of H^ decays to SM states do depend on 
a and /3, so the a • Br(gg — > A(H^ — > W^A)) ultimately varies with sin a and tan/3 due to 
the changing total width. As is apparent in Figure 11, the cross section for these processes 
is quite low, on the order of a few tens of femtobarns before further branching fractions are 
applied, so their inclusion is essentially for the sake of completeness; they contribute very 
little to the total multi-lepton limit. 

Consequently, most qualitative features of this benchmark spectrum may be understood 
simply by the combination of the direct multi-lepton decays of H and h as well as the cascade 
decay H — > AA with A — > tt, which in this spectrum is the only source of multi-lepton signals 
from processes involving the pseudoscalar. 

Types I & III 

In a Type I 2HDM, the limit is largely governed by the direct multi-lepton decays of h and H. 
In particular, the multi-lepton decays of h are SM-like around the alignment limit and decrease 
slowly away from this limit. As sin a — > —1, the multi-lepton signals of H become important 
and somewhat compensate for the vanishing signals of h. The branching ratio H —> A A is 
large at moderate sin a and large tan/3, but Br(A —> tt) does not grow exceptionally large 
in this regime, so the contribution to multi-lepton limits from H — > AA is not great. 

In the Type III 2HDM, the multi-lepton signals are much as in the Type I 2HDM with 
the exception of those involving h — > tt and A — > tt. Thus, the process gg —> H — > AA — > 4r 
contributes significantly in this 2HDM type. Unsurprisingly, in the region excluded by 5 fb _1 
data, a ■ Br(gg — )■ H — )■ AA — > 4t) is large, > 500 fb, with the current exclusion contour 
tracking the contours of T(H — > AA). 

Types II & IV 

In Type II, the multi-lepton signals of h from decays to vectors decrease rapidly above the 
alignment limit and increase rapidly below it, again supplemented by the multi-lepton signals 
of H as sin a — > —1. The multi-lepton signals of associated production with h —> tt are 
somewhat important at large tan /3, but are not significantly enhanced over the SM rate since 
h — > bb grows equally quickly and controls the total width. Similarly, although the Att 
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Figure 10. Multi-lcpton limits from the CMS multi-lcpton search with 5 fb _1 of 7 TeV proton-proton 
collisions [9] for the production and decay topologies of Benchmark Spectrum 4 given in Table 5, for 
Type I (top left), Type II (top right), Type III (bottom left), and Type IV (bottom right) couplings 
as a function of sin a and tan/3. Limits were obtained from an exclusive combination of the observed 
and expected number of events in all the multi-lepton channels presented in Table 10. The solid and 
dashed lines correspond to the observed and expected 95% CL limits on the production cross section 
times branching ratio in multiples of the theory cross section times branching ratio for the benchmark 
spectrum and 2HDM type. The blue shaded regions denote excluded parameter space. The solid 
red line denotes the alignment limit sin(/3 — a) = 1. The gray shaded region corresponds to areas of 
parameter space where vector decays of the heavy CP-even Higgs, H — > VV* , are excluded at 95% 
CL by the SM Higgs searches at 7 TeV [1]. 

coupling grows with tan f3, so too does the coupling Abb, so H — > AA — > 4r is not particularly 
important here. 
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Figure 11. 2HDM Benchmark Spectrum 4 partial width r(i/ ± — > W ± h) in units of GeV, and cross 
section times branching ratios a ■ Bi(qq —> A{H ± — > Wh)) and <t • Br(qq — > ^(i?* — > WM.)) in units 
of pb for Type I couplings. The partial width T(H ± — > W A) is independent of a and f3 and is not 
shown explicitly. 

For Type IV 2HDM the limits are much as in the Type II 2HDM, albeit with the loss of 
multi-lepton signals coming from h — > tt and A —> tt at large tan /3, leading to the weakest 
overall limits among 2HDM types. 

6 Conclusion 

In the wake of the discovery of a Standard Model-like Higgs, exploring and bounding exten- 
sions of the EWSB sector takes on paramount importance. Models with two Higgs doublets 
are among the simplest and best motivated such extensions to the Higgs sector. In this 
work, we have examined the reach of multi-lepton searches for probing the collective leptonic 
signatures resulting from the additional Higgs bosons in 2HDMs. In a study of 20 exclusive 
multi-lepton channels in four benchmark spectra with four discrete types of fermion couplings 
across 222 production and decay topologies, using a factorized mapping procedure [12] we 
determined regions of 2HDM parameter space probed by data from a recent CMS multi- 
lepton search [9] with 5 fb -1 of 7 TeV proton-proton collisions. These results provide new 
limits in some regions of 2HDM parameter space that have not been covered by other types 
of direct experimental investigations. Increased luminosity and production rates with 8 TeV 
proton-proton collisions and beyond will extend the 2HDM limits and discovery potential of 
multi-lepton searches. 

Although the CMS multi-lepton searches [8, 9] in their current incarnation are extremely 
powerful tools for probing new physics, with appropriate modifications the searches could be 
tailored in order to enhance sensitivity to 2HDM signals. Subdividing all exclusive multi- 
lepton channels by zero, one, or two or more 6-tagged jets in an event should significantly 
increase sensitivity to 2HDM final states with bottom quarks. Although many of 3- and 4- 
lepton events coming from production and decays of scalars in 2HDM populate the exclusive 
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channels with relatively high backgrounds, most of the irreducible prompt background does 
not contain additional 6-jets. For those backgrounds that do, very rarely, 6-jets will provide 
isolated leptons, so two 6-tags will substantially reduce major backgrounds (with the notable 
exception of ti plus a prompt fake lepton and ttV), while leaving many 2HDM signal processes, 
such as H ->• hh ->• ZZbb, ttA ->• tiZh, ttA ->• ttrr, H ->• ->• Z/i) ->• rrZM>, H ->• 
H + H~^ tbWh, H — >■ — > ZZ/s — >■ and, of course, it/i, relatively unaffected. 

Final states with multiple r-leptons are among the most promising for discovery or ex- 
clusion of various 2HDM. In our study, we have focused solely on leptonically-decaying rs, 
since final states with hadronic rs will often have larger backgrounds. However, ignoring 
hadronic rs reduces sensitivity to, in particular, four-r final states with low a ■ Br. A further 
partitioning of the 4£, 2r bins in a study optimized for four-r signals may yield lower back- 
grounds in DYO bins, e.g. r^r^e - ^ - , allowing for improved limits. As much of the energy 
in these events are going into leptons, defining signal regions either with harder pt cuts on 
leptons or with a cut on ^ p^/ could serve to significantly deplete the high SM backgrounds 
in some bins while leaving the signal largely unfazed. We have also restricted our focus to 
three- and four-lepton final states. Some additional sensitivity may be gained by adding 
exclusive channels with same-sign di-leptons subdivided by various combinations of $ T and 
Ht- These channels would capture other decay modes of some of the production and decay 
topologies studied here, as well as bring in additional topologies that do not yield three or 
more leptons. Multiple Higgs bosons can also give rise to rare five- or more lepton signatures; 
adding channels to separate out these signatures would also increase sensitivity, particularly 
at high luminosity. 

Finally, with a known Higgs mass, one can capitalize on partial or full kinematic con- 
straints of its decays to help to isolate Higgs particles arising via new sources of associated 
production. Such kinematic tagging can serve to further reduce SM backgrounds. One exam- 
ple of this would be forward jet tagging to highlight VBF signals. Another would be channel 
specific lepton kinematics focussed at specific decay topologies. One of the simplest and most 
effective ways to utilize kinematic tagging to enhance sensitivity to certain multi-lepton signa- 
tures that include a SM-like Higgs boson would be to subdivide the DY2 four- or more lepton 
channels into an On Higgs category in which the invariant mass of the four leptons fall within 
a small window centered on the Higgs boson mass. Signals that include at least one SM-like 
Higgs boson that decays directly to four leptons fall in this sub-channel. The backgrounds 
in this special On Higgs sub-channel are very limited, thereby increasing sensitivity to such 
Higgs boson signals. Utilizing partial (rather than full) kinematic tagging could also increase 
sensitivity to other decay topologies that fall in other channels. 

While we have focused on 2HDMs, other extensions of the Higgs sector can lead to 
the production of new heavy, Higgs-like scalar resonances with decay topologies similar to 
those studied in this work. Such new, Higgs-like particles generally lead to intermediate 
states composed of the heaviest SM particles, including t, h, Z, W, b and r, whose final 
states contain multi-lepton signatures. If there exists an extended Higgs sector, multi-lepton 
searches optimized for the leptonic final states of Higgs scalars may prove an effective route 
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for discovering new physics beyond the Standard Model. 
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